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ABSTRACT 

The  application  of  bonded  repairs  to  structures  subject  to  acoustic  fatigue  is 
investigated,  with  a  view  to  developing  simple  solutions  suitable  for  engineering 
applications.  Acoustic  fatigue  is  caused  by  a  random  response  of  a  structure  due  to 
time  varying  pressure  waves  from  engine  and/or  aerodynamic  sources.  For  vibration 
problems  the  use  of  damping  in  repairs  will  reduce  the  amplitude  of  vibration  and 
hence  the  stress  intensity  factor  in  the  cracked  structure.  To  assess  the  viability  of  the 
composite  bonded  repairs  subject  to  acoustic  fatigue  new  analytical  tools  are  needed. 
In  this  paper  the  use  of  a  simplified  dynamic  analysis  combined  with  crack-bridging 
theory  allows  the  influence  of  damping  to  be  assessed.  This  theory  is  adapted  to 
panels  with  various  edge  support  conditions  with  good  results.  Validation  of  this  work 
is  carried  out  using  a  Finite  Element  (F.E.)  analysis  in  which  three-dimensional  brick 
elements  are  used  to  model  the  plate,  adhesive  and  boron  patch,  and  use  is  made  of  the 
power  spectral  density  capability  of  the  NASTRAN  program. 
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Computation  of  the  Stress  Intensity  Factor  in 
Bonded  Repairs  to  Acoustically  Fatigued  Panels 


Executive  Summary 

The  repair  of  cracked  aircraft  structures  subject  to  manoeuvre-type  loads  using 
composite  bonded  repairs  has  resulted  in  considerable  aircraft  life  time  extension  and 
hence  cost  savings.  However  the  use  of  bonded  patches  to  repair  panels  with 
acoustically  induced  cracks  (acoustic  fatigue)  is  only  recent  and  presents  some 
difficulties.  Acoustic  fatigue  is  a  result  of  high  frequency  lateral  vibration  of  an  aircraft 
panel  as  a  result  of  intense  random  noise  caused  by  engine  and/or  aerodynamic 
effects.  For  example  acoustically  induced  cracks  have  been  found  in  the  lower  external 
surface  of  the  nacelle  skin  of  the  F/A-18  aircraft.  In  this  case  the  overall  sound  pressure 
levels  of  the  order  of  170  db  have  been  recorded.  Attempts  to  repair  these  cracks  by 
applying  standard  designs  of  composite  bonded  repair  were  made,  however  the  cracks 
continued  to  grow.  These  panels  were  repaired  on  the  basis  of  high  level  manoeuvre 
loads  but  low  frequency  out-out-of  plane  loads.  It  is  evident  that  the  use  of  composite 
bonded  repairs  to  cracked  panels  subject  to  acoustic  excitation  requires  analytical  tools 
that  take  high  frequency  random  out-of-plane  vibration  into  account. 

Previous  work  has  shown  that  while  the  boron  fibre  composite  bonded  repairs  applied 
to  panels  containing  acoustically  induced  cracks  do  reduce  the  driving  force  at  the 
crack  (stress  intensity),  the  high  number  of  cycles,  due  to  the  high  loading  frequencies, 
lead  to  significant  crack  growth.  For  example,  one  hour  of  flying  may  result  in  about 
106  cycles  for  the  high  frequency  acoustic  fatigue  case.  Since  it  is  known  that  the 
amplitude  of  vibration  of  the  panel  is  dependent  on  the  damping,  then  increasing  the 
damping  may  reduce  the  crack  growth  rate  to  an  acceptable  level. 

In  this  report  an  analytical  solution  is  proposed  to  allow  the  computation  of  stress 
intensity  in  highly  damped  repairs,  which  involves  the  use  of  crack-bridging  theory. 
This  analytical  approach  is  validated  against  detailed  3D  finite  element  analysis. 
Furthermore  a  simplified  dynamic  analysis  of  the  random  loading  is  adopted.  It  is 
shown  that  a  combination  of  the  simplified  dynamic  analysis  and  crack  bridging 
theory  is  applicable  to  infinitely  wide  panels.  It  has  been  shown  in  previous  work  that 
a  simple  plate  analysis  for  crack  bridging  can  predict  the  behaviour  for  plates  in 
cylindrical  bending  for  the  static  case.  Furthermore  use  of  simple  plate  bending  theory 
can  extend  this  theory  to  rectangular  panels  restrained  on  all  edges.  As  a  result  this 
theory,  together  with  appropriate  crack  growth  data  and  damping  data,  will  enable  the 
design  of  highly  damped  repairs  of  acoustically-induced  cracked  panels  of  aircraft 
structures. 
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1  Introduction 


The  repair  of  cracked  aircraft  structures  subject  to  in-plane  loads  using  bonded  repairs 
has  resulted  in  considerable  aircraft  life  time  extension  and  hence  cost  savings  [1]. 
However  the  use  of  bonded  patches  to  repair  panels  with  acoustically  induced  cracks 
(acoustic  fatigue)  is  only  recent  and  presents  some  difficulties.  Acoustic  fatigue  is  a 
result  of  high  frequency  lateral  vibration  of  an  aircraft  panel  due  to  time  varying 
pressure  waves  caused  by  engine  and/or  aerodynamic  effects.  For  instance, 
acoustically  induced  cracks  have  been  recorded  in  the  lower  external  surface  of  the 
nacelle  skin  of  the  F/A-18  aircraft,  as  shown  in  Fig.  1.  In  this  case  the  overall  sound 
pressure  levels  of  the  order  of  170  db  have  been  recorded,  [2]  and  [3].  Attempts  to 
repair  these  cracks  by  applying  standard  methods  of  bonded  repair  [4]  developed  for 
in-plane  loads  were  made,  however  the  cracks  continued  to  grow.  The  application  of 
bonded  repairs  to  acoustically-induced  cracks  requires  analytical  tools  that  take  into 
account  high  frequency  out-of-plane  vibration.  It  has  been  found  [5]  that  while  the 
boron  fibre  composite  bonded  repairs  did  reduce  the  stress  intensity,  at  the  crack-tip, 
the  high  number  of  cycles  lead  to  significant  crack  growth  of  up  to  2 mm  per  flight  hour. 

Since  it  is  known  that  the  amplitude  of  vibration  is  inversely  proportional  to  the  square 
root  of  the  damping,  then  a  combination  of  damping  and  stiffness  may  reduce  the 
crack  growth  rate  significantly.  The  application  of  highly  damped  repairs  (Durability 
Patches)  to  acoustically  damaged  panels  have  been  proposed  by  [6].  Furthermore  [7] 
have  applied  a  highly  damped  repair  to  the  vertical  fin  of  the  F-15  aircraft,  as  part  of  a 
test  program. 

In  this  paper  an  analytical  solution  will  be  developed  to  allow  the  computation  of 
stress  intensity  factors  in  highly  damped  repairs,  which  involves  the  use  of  crack- 
bridging  theory  [8].  A  correction  factor  has  been  introduced  to  extend  the  crack 
bridging  theory  to  finite  width  rectangular  panels. 


Figure  1:  Location  of  the  cracking  in  the  lower  nacelle  inlet. 
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2  Analytical  Approach 

Acoustic  fatigue  is  associated  with  the  random  loading  of  a  structure  over  a  wide 
frequency  range.  In  the  case  considered  here  the  loading  is  applied  to  a  rectangular 
panel  shown  in  Fig.  2,  with  three  different  sets  of  boundary  conditions.  The  complete 
analytical  solution  of  the  dynamic  response  of  the  panel  is  very  complex  [9].  However 
in  the  case  of  acoustic  fatigue  a  number  of  simplifying  assumptions  can  be  made  as  the 
emphasis  here  is  to  determine  the  crack  growth  driving  force.  Firstly  the  time  varying 
pressure  can  be  considered  to  be  uniform  and  in  phase  over  the  panel.  Secondly,  since 
the  first  mode  contributes  to  98%  of  the  panel  response,  [3],  then  the  contribution  of  all 
other  modes  will  be  ignored.  Also  the  variation  of  power  spectral  density  (FSH)  may  be 
assumed  to  be  constant  with  frequency  near  the  first  mode  frequency.  Based  on  these 
assumptions  the  root-mean-square  of  the  stress  response  is  given  [10,11]  as: 


cr 


b.rms 


tfoSAfo) 


4  8 


<y , 


(l) 


where  <Jb  rms  is  the  resultant  root-mean-square  stress  at  any  point  &  in  the  plate  (MPa) 
crb  is  the  stress  at  point  b  due  to  a  static  unit  pressure  (MPa/Pa) 

Sj  ( f0  )  is  the  FSD  of  the  excitation  at  f0  (Pa 2  /  Hz) 

fQ  is  the  first  resonant  frequency  (Hz)  which  depends  on  boundary 
conditions 

8  is  the  viscous  damping  ratio  for  the  overall  structure 

A  similar  expression  for  equation  1  [9]  can  be  derived  for  displacements. 

Since  the  stress  intensity  is  proportional  to  the  prospective  stress  (<7rmi)  prevailing  in 
the  uncracked  plate  subjected  to  random  acoustic  pressure,  then  the  following  equation 
also  holds: 


Kms  = 


¥oS,(fo) 


48 


K 


(2) 


where  Krms  is  the  root-mean-square  (in  time  space)  of  the  stress  intensity  factor  and  K 
is  the  stress  intensity  factor  due  to  the  application  of  a  unit  pressure  on  a  cracked 
structure  with  or  without  repairs.  For  a  repaired  structure,  K  will  be  evaluated  using 
crack-bridging  theory.  The  problem  is  idealised  as  shown  in  Fig.  2a  in  which  the  plate 
is  simply  supported  along  two  edges  (cylindrical  bending)  subject  to  a  random  loading. 
The  first  resonant  frequency  is  given  by: 
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1C  D 


where  D  = 


12(1  -o2) 

E  is  Young's  modulus, 
tp  is  the  thickness  of  the  plate, 
p  is  the  density , 

L  and  b  are  the  length  and  width  of  the  plate  respectively. 


For  the  case  of  a  plate  with  all  edges  restrained,  Fig.  2b  and  2c,  the  first  resonant 
frequency  is  given  by: 


/o=V-/(2*/£)J-^7 

2n  \  pt  L 


and  f(2b/  L)  is  given  in  [12]  for  either  or  all  simply  supported  or  all  clamped  edges 
and  is  shown  in  Fig.  3. 


(a)  simply  supported  on  2  edges  (b)  simply  supported  on  all  edges 


(  c  )  clamped  on  all  edges 


Figure  2;  Edge  support  conditions  for  panels. 
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Figure  3:  Frequency  function  for  first  mode  of  a  rectangular  plate,  all  edges  simply  supported 
or  all  edges  clamped. 


While  in  practise  most  acoustic  fatigue  cracking  occurs  at  the  edge  of  panels,  for 
simplicity  the  crack  has  been  assumed  to  be  at  the  centre  of  the  panel  as  shown  in 
Fig.  2.  For  the  case  of  simply  supported  edges,  the  bending  moment  that  is  experienced 
by  the  crack  at  the  centre  of  the  plate  is  higher  than  that  at  the  edge.  This  location  will 
also  simplify  the  use  of  crack  bridging  theory.  The  size  of  the  panel  consided  here  is 
b=126mm  and  L=195mm. 


3  Sound  Pressure  Levels 

In  this  case,  the  sound  pressure  levels  (SPL)  measured  on  the  external  surface  of  F/  A- 
18  inlet  nacelle  will  be  used,  and  are  plotted  in  Fig.  4  with  increasing  frequency.  The 
spectrum  level,  relative  to  the  overall  sound  pressure  level  (OASPL),  was  derived  from 
in-flight  one  third  octave  SPL  measurements,  [2].  This  spectrum  level  is  now  used  to 
calculate  S,  (/0 ) .  The  relationship  between  the  spectrum  SPL  and  the  r.m.s.  fluctuating 
pressure  p,  (units  are  Pa),  is  given  by  [13]  as: 

Prms=  io(S/>i/20-4-7)  (5) 


4 


DSTO-TR-1061 


and  the  power  spectral  density  of  acoustic  pressure,  i.e.  PSD  of  the  excitation,  at  any 
given  frequency  is  given  by: 


10^™  (6) 

The  curve  in  Fig.  4  has  been  approximated  by  two  straight  lines  defined  by  the  three 
points  listed  in  Table  1  and  shown  in  Fig.  4.  This  approximate  spectrum  is  used  as  the 
excitation  PSD  for  the  F.E.  analysis. 


Table  1:  Input  power  spectral  density. 


Frequency  (Hz) 

Three  points  to 
appromimate 
curve  (dB) 

Pressure  spectrum 

level  (dB) 

Wo) 

(MPa)2/Hz 

31.5 

-32.2 

140 

4.0  xlO-8 

1000. 

-35.2 

137 

2.005x10-8 

8000. 

-48.1 

124.1 

1.028xl0-9 

Frequency  (Hz) 

Figure  4:  Spectrum  and  one-third  octave  band  levels  of  sound  pressure  over  the  external 
nacelle  inlet  (where  OASPL=170  db)  [2]. 
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4  Crack  Bridging  Theory 

The  most  common  type  of  repairs  to  be  applied  to  aircraft  skins  are  one  sided  as 
illustrated  in  Fig.  5  since,  in  general,  access  is  only  available  to  the  outside  of  the  skin. 
Analytical  solutions  for  one  sided  repairs  are  available  for  both  in-plane  [14]  and  out- 
of-plane  loading  [8].  A  significant  feature  of  bonded  repairs  is  that  the  variation  of 
stress  intensity  K(z )  with  crack  length  does  not  increase  indefinitely  with  crack 
length.  Instead,  the  stress  intensity  approaches  an  asymptotic  value.  The  thickness  of 
the  adhesive  is  included  in  the  analysis  to  represent  the  damping  material.  Crack 
bridging  analysis  for  the  case  shown  in  Fig.  5b  gives  the  solution,  shown  in  Fig.  6,  as: 

K(z)  =  ab.fm+ab~fb  (7) 

lp 

where  fm  and  fb  denote  the  membrane  and  bending  components  of  the 
normalised  stress  intensity  factor  induced  by  a  uniform  bending  stress,  which 
varies  linearly  through  the  plate  thickness  [15].  ab  is  the  maximum  bending 
stress  in  absence  of  a  repair  and  crack.  Both  fm  and  fb  depend  on: 

Ep ,  Young's  modulus  for  the  skin, 

Er  ,  Young's  modulus  for  the  reinforcement, 
tp ,  the  thickness  of  the  skin, 

ta ,  the  thickness  of  the  adhesive  or  damping  material, 
tr ,  the  thickness  of  the  reinforcement, 

L,  the  length  of  the  beam. 
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(a)  Un-cracked  case 


h - l  - H 


(b)  Cracked  and  repaired  case 


Figure  5:  One  sided  repair  consisting  of  a  reinforced  beam  subject  to  one  unit  of  lateral 
pressure. 


Table  2:  Material  properties. 


Material 

Young's 
modulus  (GPa) 

Poisson's  ratio 

Density 
(Mg/  mm3) 

Aluminium 

71. 

.33 

2.77x10-9 

Adhesive  [16] 

2.3 

.35 

1.2x10-9 

Boron  [17] 

207. 

.21 

2.0x10-9 

The  crack-bridging  model  involves  a  reinforced  plate  with  a  central  crack,  as  shown  in 
Fig.  5b.  It  is  worth  noting  that,  for  a  unit  pressure,  ub  needs  to  be  calculated  using  the 
appropriate  boundary  conditions,  e.g.  simply  supported  or  clamped.  The  crack¬ 
bridging  model  plots  of  K{z)l  ob  in  Fig.  6  show  the  maximum,  minimum  and  mean 
values  of  stress  intensity,  ( Kmax/  Kmean ,  and  Kmtn).  In  this  case,  the  material  properties  and 
dimensions  are  given  in  Table  2.  The  stress  cxb  is  the  bending  stress  at  the  centre  of  the 
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beam,  in  the  absence  of  a  repair  and  a  crack.  Note  that  Fig.  6  shows  the  expected 
asymptotic  behaviour. 


Figure  6:  Results  from  crack  bridging  model,  [Wang  and  Rose]  for  constant  lateral  static 
loading  of  a  beam. 


Firstly  the  mean  stress  intensity  factor  is  given  by: 
K  4 .  k 

jy-  _  ^  max _ min 

^  mean  2 

and  the  bending  stress  intensity  factor  is  given  by: 
K  —  K 

zr  ^  max  1Vmin 

Kb~  2 

These  equations  will  be  considered  futher  on. 


(8) 


(9) 


For  a  simply  supported  plate.  Fig.  2a,  subject  to  a  uniform  pressure  of  unit  magnitude, 
the  maximum  bending  stress  is  given  by: 
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cr 


b 


3  L2 


(10) 


In  the  context  of  repairs  ob  corresponds  to  the  nominal  far-field  bending  stress.  Hence 
using  equations  (1  and  10)  the  root-mean-square  bending  stress  of  a  simply  supported 
plate  along  two  edges  due  to  the  time  varying  random  pressure  with  PSD,  Sf  ( f  )  is 
given  by: 


b.rms 


|joW  3L 

4 5  4 tl 


(11) 


5  Finite  Element  Approach 


The  random  response  analysis  capability  of  the  NASTRAN  program  [18]  has  been  used 
to  validate  the  analytical  results  presented  here.  The  system  response  to  random 
excitation  is  evaluated  via  the  frequency  response  technique  viz., 

£,(<»)  (12) 

where  in  this  case  H(co )  is  the  frequency  response  of  a  displacement  us  due  to  an 
excitation  force  (pressure)  pb .  Sj  and  S;  are  the  PSD  of  the  response  (displacement) 
and  of  the  excitation  (pressure),  respectively.  This  analysis  also  allows  the  statistical 
properties  of  the  system  to  be  evaluated.  Random  vibrations  considered  here  involve 
all  frequencies  at  any  one  instant  in  time.  After  calculating  the  PSD,  the  root-mean- 
square  (r.m.s.)  of  the  response  can  be  determined  by  computing  the  square  root  of  the 
PSD  area: 


Jrms  d0)  (13) 

where  the  limits  of  the  integral  will  be  truncated  as  shown  in  Fig.  4.  No  structural 
response  is  expected  below  30Hz  and  the  excitation  levels  above  8000Hz  are  minimal. 

A  similar  application  of  finite  element  techniques  to  undertake  a  PSD  analysis  to 
acoustic  fatigue  problems  has  been  reviewed  by  [13]. 
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5.1  Computation  of  stress  intensity  factors  using  F.E. 


In  the  F.E.  model  the  thickness  of  the  plate  is  modelled  using  three  layers  of  20  noded 
brick  elements  in  order  to  account  for  the  bending  behaviour  of  the  skin.  The  skin 
thickness  is  approximately  1  mm,  hence  the  condition  of  plane  stress  is  assumed.  The 
computation  of  the  stress  intensity  factor  is  determined  directly  from  the  crack  opening 
displacements  near  the  crack  tip  using  displacements  from  the  PSD  analysis.  The  r.m.s. 
crack  tip  stress  intensity  factors  for  mode  l  are  derived  from  the  standard  asymptotic 
relation  [20]: 


where  Krm  is  the  root-mean-square  stress  intensity  factor  evaluated  at  the  desired 
through  thickness  location 

E  is  Young's  modulus 

Urms  is  the  root-mean-square  of  the  crack  opening  displacement  at  distance  l 
from  the  crack  tip  (perpendicular  to  the  plane  of  the  crack)  and  at  the  desired 
through  thickness  location. 

The  F.E.  model  is  comprised  entirely  of  20  noded  brick  elements.  For  the  plate  three 
elements  are  used  through  the  thickness  while  two  elements  are  used  through  the 
thickness  for  the  adhesive  and  patch.  The  size  of  the  crack  tip  element  (/  =  0.15625 mm ) 
has  been  selected  such  that  /  «  1  /  k ,  where  k  =  (crb/  ) 2  »  1 .0 mm ,  where  values  of 

K/ab  are  taken  from  Fig.  6. 


6  Analysis  of  Results 

6.1  Plate  with  two  edges  simply  supported 

6.1.1  Plate  with  no  crack  and  no  repair 

Consider  the  configuration  shown  in  Fig.  2a  in  which  the  plate  is  simply  supported  at 
two  edges  but  does  not  contain  a  crack  or  a  repair.  The  material  properties  of  this  plate 
are  shown  in  Table  2  while  the  dimensions  are  shown  in  Table  3.  As  an  approximation 
the  stress  crb  at  x=0  of  the  plate  subject  to  unit  pressure  is  given  by  equation(lO).  The 
frequency  of  die  first  mode  (69  Hz)  is  computed  using  equation(3).  The  corresponding 
PSD  of  the  excitation  at  f0,  Sb (/„)  is  interpolated  from  the  spectrum  listed  in  Table  1, 
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and  is  shown  in  Table  4.  Using  equation(l)  the  rms  stress  ( crrms )  is  evaluated,  with  a 
viscous  damping  ratio  of  0.108  at  the  centre  and  on  the  surface  of  the  plate  and  is  given 
in  Table  5.  Also  listed  in  Table  5  is  the  F.E.  result  which  in  this  case,  shows  a 
disagreement  of  4.4%  compared  to  the  analytical  bending  stress  results. 


Table  3:  Dimensions  of  plate  shown  in  Fig.2. 


L  (mm) 

2b  (mm) 

tp  (mm) 

ta  (mm) 

tr  (mm) 

195. 

252 

1.143 

0.254 

0.635 

Table  4:  Dynamic  properties  of  cracked  and  repaired  plates 


No.  edges  restrained 

fo  (Hz) 

Si  (f0 )  (MPa)2/ Hz 

5 

2 

69  (equation  3) 

3.917X10-8 

.108 

4 

118  (equation  4) 

3.822x10-8 

.108 

Table  5:  Stress  results  for  un-cracked,  un-repaired  plate. 


Theoretical  (MPa) 

F.E.  result  (MPa) 

Error  % 

96.8 

92.7 

4.4 

6.1 .2  Plate  with  crack  and  repair 

In  this  case  the  properties  of  the  repaired  plate  are  also  contained  within  Tables  2 
and  3.  The  patch  is  assumed  to  cover  the  full  width  of  the  plate.  The  following 
equation,  derived  from  equations(2  and  8),  has  been  used  to  calculate  the  root-mean- 
stresses  where  values  of  /  ab ,  /  oh  are  given  in  Fig.  6: 


RW„)  3iJ 

iK) 

\  45  4 1) 

(15) 


where  the  frequency  is  given  by  equation(3). 

Using  equations(8  and  9)  the  ratio  of  {Kmem)rm  l{Kb)rms  has  been  computed  for  both 
F.E.  and  analytical  analysis  and  is  shown  in  Fig.  7.  The  best  agreement  occurs  at  a  half 
crack-length  of  approximately  20  mm  and  corresponds  to  an  error  of  6.7%.  The  largest 
error  occurs  at  the  smallest  and  largest  crack  length.  An  average  error  between  these 
results  is  9%.  The  characteristic  of  the  F.E.  results  is  the  constant  value  of 


11 


DSTO-TR-1061 


(K mean ) rms  l(.Kb)rms  as  crack  length  increases  while  in  the  case  of  theory  the  ratio  has 
yet  to  stabilise.  Also  the  closed  form  solution  results  are  always  higher  than  the  F.E. 
results. 


a  (mm) 


Figure:  7:  Comparison  of  (Kmean )rms  /(Kb)ms  for  theory  and  F.E. 


0.0  0.2  0.4  0.6  0.8  1.0 


Figure  8:  Plot  of  (Kb)rmsF.E./(Kb  )rms  Theory ,  polynomial  fit:  l.+0.011(a/b)+0.366(a/b)2 

This  error  may  be  due  to  an  edge  effect  as  a/b->  1.  In  order  to  make  a  correction  factor 
(Kb)rmf  /(Kb)lr  has  been  plotted  against  a/b,  as  shown  in  Fig.  8,  and  fitted  with  a 
second  order  polynomial  such  that  the  results  are  assumed  to  be  correct  at  small  crack- 


12 


DSTO-TR-1061 


lengths.  It  is  noted  that  the  ratio  (Kmean)Fmf  /(K mean)™esory  should  always  be  one.  This 
polynomial  is  defined  below  where  the  values  of  a  and  b  are  defined  in  Fig.  2: 

/  =  (Jf»  )"■  KKt  )2r  =  1.0+0.011(a/b)+0.366(a/b)2  (16) 


As  a  result  the  following  correction  is  applied: 

( V-  \  _  \Theory 

K^b/rms  ~  J  \*-b)rms 

(17) 

(18) 

(19) 

where  )£f  =  (Km„, 

A  comparison  with  the  corrected  theory  and  F.E.  results  for  )rms ,  (Kmean  )rms  and 
(Kmn )rms  are  shown  in  Fig.  9.  The  corrected  values  of  (K^)rms,(Kmean)rms  and 
(K^  )rms  for  crack  bridging  theory  are  all  conservative. 


Figure  9:  Comparison  of  analytical  and  F.E.  results  for  two  edges  simply  supported. 
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In  the  analytical  solution  it  is  assumed  that  the  frequency  does  not  change  with  crack 
length.  The  F.E.  values  of  first  resonant  frequency  with  crack  length  ratio  are  shown  in 
Table  6,  and  indicate  that  only  a  small  change  of  frequency  occurs,  (i.e.,  frequency 
change  of  1  %  compared  to  the  original  first  resonant  frequency  for  crack  length  ratio  of 
a/b=.92). 

Table  6  Frequency  of  repaired  plate  versus  crack  length  ratio  for  two  edges  simply  supported. 


Crack  length  ratio  a/b 

Frequency  Hz 

.0397 

72.43 

.0794 

72.41 

.1587 

72.33 

.3175 

72.13 

.4762 

71.90 

.6349 

71.64 

7937 

71.36 

.9127 

71.14 

6.2  Plates  simply  supported  on  all  edges 

In  general  aircraft  skin  panels  are  usually  connected  to  other  plates  or  stiffeners  at  their 
boundaries,  consequently  the  acoustic  response  lies  somewhere  between  the  limiting 
cases  of  simply-supported  edges  or  fully-fixed  edges.  Simple  plate  theory  can  be  used 
to  determine  the  bending  moment  distribution  in  the  plate  resulting  from  a  uniform 
pressure  distribution.  In  the  case  of  four  simply  supported  edges  the  bending 
distribution  relative  to  the  crack  location  is  shown  in  Fig.  10,  where  the  plate  is  of 
dimensions  2b  by  L.  The  average  bending  moment  can  be  determined  from  [20]  as: 

Mme=pL2  (20) 


Note  that  /?  is  a  function  ot{a/b  and  b/L)  and  is  shown  in  Fig.  11. 

Since:  <Jb  =  6M/t2  then  from  equation(20): 

ab=6/3L2/t2p  (21) 

The  stress  intensity  factor  is  approximately  equal  to  that  corresponding  to  a  uniform 
bending  moment  with  magnitude  being  equal  to  the  average  over  the  crack  length. 
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Hence  from  equation(2): 


Wo  Wo)  6/£2 

i  4^  t] 

(22) 


where  K(z )  is  given  by  equation(7).  In  this  case  the  frequency  f0  is  given  by 
equation(4)  and  S/fo)  is  interpolated  from  Table  1,  and  is  shown  in  Table  5. 


Figure  10  Comparison  of  bending  moments  from  beam  and  plate  with  simply  supported  edges. 

The  closed  form  result  of  ( K mx)rms,  shown  in  Fig.  12,  gives  conservative  agreement 
over  the  complete  crack  length  range.  It  is  evident  in  the  problem  considered  that  for 
an  aspect  ratio  of  2&/L«  1,  the  edges  (whose  dimension  is  L  as  shown  in  Fig.  10)  have  a 
strong  restraining  influence  on  the  bending  developed  in  the  repair,  i.e.  the  F.E.  results 
show  lower  values  than  those  from  the  closed  form  solution.  For  large  aspect  ratio's 
(2b/L»l)  it  would  be  expected  that  closer  agreement  would  occur.  A  closed  form 
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6.3  Plates  clamped  on  all  edges 

The  F.E.  results  for  built  in  edges  are  also  shown  in  Fig.  12  however  an  appropriate 
correction  factor  has  not  been  applied  to  crack  bridging  theory.  The  solution  is  also 
applicable  for  cracks  along  the  boundary,  which  generally  occur  in  practice,  however 
the  result  will  be  conservative  since  the  bending  moment  will  be  lower  at  the  boundary 
than  that  at  the  centre  of  the  plate. 

6.4  Effect  of  damping 

For  the  case  of  a  plate  with  all  edges  simply  supported  the  closed  form  solution 
has  been  used  to  consider  the  influence  of  damping.  Results  shown  in  Fig.  13,  are  a 
plot  of  with  viscous  damping  coefficient,  and  demonstrates  the  possible 

reduction  in  stress  intensity  factor  with  increased  damping.  The  low  values  of 
damping  would  correspond  to  a  metallic  riveted  aircraft  skin  [2],  and  the  high  values 
to  a  panel  treated  with  high  damping.  It  is  evident  that  for  a  highly  damped  repair  the 
stress  intensity  factor  may  be  further  reduced  by  a  factor  of  2.5. 


Figure  13:  Variation  of  Stress  intensity  factor  with  increasing  viscous  damping  coefficient. 
(Note  the  reduction  of  stress  intensity  factor  for  highly  damped  repair  in  comparison  to  metallic 
riveted  structure). 

6.5  Power  spectral  density  of  the  response  (PSD) 

The  PSD  of  the  response  using  F.E.  is  shown  in  Fig.  14  for  the  case  of  a  highly  damped 
repaired  plate  simply  supported  on  all  edges.  Since  the  response  is  given  by 
equation(12),  then  a  running  integral  of  the  PSD  of  the  response  will  give  the 
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contribution  of  each  mode.  A  COD  method  is  used  to  determine  the  stress  intensity 
using  displacements  at  the  crack  tip.  It  can  be  seen  that  the  first  mode,  at  a  frequency  of 
118  Hz,  dominates  the  response  of  the  structure  compared  to  the  higher  modes.  A 
running  integral  of  this  curve  shows  that  the  first  mode  contributes  to  95%  of  the 
response.  While  there  are  extreme  cases  in  which  this  does  not  occur  [13],  in  general 
however,  the  first  mode  does  contribute  to  most  of  the  response,  and  hence  this 
assumption  made  to  derive  the  analytical  solutions  for  is  accurate  to  the  degree 
shown  in  this  work. 
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Figure  14:  The  PSD  of  the  displacement  response  for  a  plate,  with  a  highly  damped  repair, 
simply  supported  on  all  edges. 


7  Conclusions 


In  this  paper  an  analytical  solution  has  been  developed  using  a  simplified  dynamic 
analysis  together  with  crack  bridging  theory  to  give  a  useful  upper  bound  for  idealised 
configurations.  It  has  been  shown  that  the  analysis  for  crack  bridging  can  predict  the 
behaviour  for  plates  in  cylindrical  bending  and  rectangular  panels  restrained  on  all 
edges.  Furthermore  the  solution  is  also  applicable  for  plates  containing  cracks  on  the 
boundary,  and  this  solution  will  be  conservative.  As  a  result  this  theory  together  with 
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appropriate  crack  growth  data,  and  damping  data  will  aid  in  the  design  of  highly- 
damped  repairs  of  acoustically-induced  cracked  panels. 
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